Effect of surface anisotropy on the hysteretic properties of a magnetic particle 
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We study the influence of surface anisotropy on the zero-temperature hysteretic properties of a 
small single-domain magnetic particle, and give an estimation of the anisotropy constant for which 
deviations from the Stoner-Wohlfarth model are observed. We consider a spherical particle with 
simple cubic crystalline structure, a uniaxial anisotropy for core spins and radial anisotropy on the 
surface, and compute the hysteresis loop by solving the local Landau-Lifshitz equations for classical 
spin vectors. We find that when the surface anisotropy constant is at least of the order of the 
exchange coupling, large deviations are observed with respect to the Stoner-Wohlfarth model in the 
hysteresis loop and thereby the limit-of-metastability curve, due to the non-uniform cluster-wise 
reversal of the magnetisation. 

I. STATEMENT OF THE PROBLEM 

In this work, we deal with the effect of surface anisotropy on the hysteretic properties (hysteresis loop and limit- 
of-metastability curve, the so-called SW astroid |gj) of a single-domain spherical particle of simple cubic crystalline 
structure, uniaxial anisotropy in the core, and radial single-site anisotropy for spins on the boundary. We compute the 
hysteresis loop and thereby the critical field of this particle by solving, at zero temperature, the local Landau-Lifshitz 
equation (LLE) derived from the classical anisotropic Dirac-Heisenberg model in field (|l|). Doing this for different 
values of the polar angle ip between the applied field and the core easy axis, renders the critical field as a function of 
ip, that is the SW astroid. Finally, our model Hamiltonian reads || 

M 

where Si is the unit spin vector on site i, H is the uniform magnetic field applied in a direction tp with respect to the 
reference z axis, AT is the total number of spins; J > is the nearest-neighbour exchange coupling; the last term in 
(|l|) is the uniaxial anisotropy energy with easy axis e* and constant Ki > 0. Core spins have an easy axis along z and 
constant K c , while surface spins have radial easy axes and constant K s . In the sequel, we use the reduced parameters, 
j = J/K c ,k s = K s /K c . 

The only free parameter of our model is the surface anisotropy constant k s . So upon varying it we estimate its 
value at which deviations from the (macroscopic) Stoner-Wohlfarth model are observed. Details of the method used 
here and more results can be found in Ref. || ■ 

II. RESULTS AND DISCUSSION 

Fig. [j] shows that when k s becomes comparable with j, the competition between exchange coupling and surface 
anisotropy produces large deviations from the SW model. Namely, the hysteresis loop exhibits multiple jumps, which 
can be attributed to the switching of different spherical shells of spins starting from surface down to the centre. The 
hysteresis loop is then characterised by two field values: One that marks the limit of metastability, called the critical 
field or the saturation field, and the other that marks the magnetization switching, and is called the switching field or 
the coercive field. In Fig. |^a we present the variation with the particle's diameter of the critical field (in diamonds) 
obtained from the numerical solution of the LLE for j = 10 2 , and (in circles) the SW critical field multiplied by the 
core-to-volume ratio, i.e. N c /N, analytically obtained upon assuming infinite exchange coupling. This figure also 
shows that for such a value of k s , when the field is applied along the core easy axis, the critical field increases with the 
particle's size, and that for k s = 1 = 10~ 2 j all hysteresis loops can be scaled with those rendered by the SW model. 
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FIG. 1. Hysteresis loop for ij) = 0, k 3 — 1 and different values of j. N = 360. 




FIG. 2. a) (in diamonds) Switching field for k s = 1, j = 10 2 versus the particle's diameter N. (in circles) SW switching 
field multiplied by N c /M . b) Astroid for k s = 1, j = 10 2 for different values of the surface-to-volume ratio N st = N s /Af. 

On the other hand, these results clearly show that even in a large particle there still exist some magnetic inhomo- 
geneities that lead to deviations from the SW hypothesis of uniform rotation, as this is only rigorously realized in an 
infinite system. Moreover, Fig. || shows that the critical field of a spherical particle with k s = 1 can be obtained from 
the SW model through a scaling with constant N c /Af. One should also note that the astroid for all particle sizes falls 
inside that of SW, see Fig. |^b. Therefore, for k s /j ~ 0.01 our model renders hysteresis loops and limit-of-metastability 
curves that scale with the SW results for all values of the angle ip, the scaling constant being N c /Af, which is smaller 
than unity for a particle of any finite size. On the other hand, the critical field increases with the particle's size and 
tends to the SW critical field in infinite systems. 



2 




FIG. 3. a) Astroid for j — 10 2 , Af = 360 and different values of surface anisotropy constant k s . The full dark line is the SW 
astroid scaled with N c /N ', but the dotted line is only a guide for the eye. b) Switching field versus the particle's diameter TV 
for V = 0, j = k s = 10 2 . 

For larger values of k s /j, but k s / j <0.2, we still have some kind of scaling but the corresponding constant now 
depends on if>. This is reflected by a deformation of the limit-of-metastability curve. More precisely, as shown by 
Fig. |3|a, the latter is depressed in the core easy direction and enhanced in the perpendicular direction. However, there 
is still only one jump in the hysteresis loop implying that the magnetization reversal can be considered as uniform. 
On the other hand, Fig. ^|b shows that for k s /j ~ 1 there are two new features in comparison with the previous 
case (compare Fig. ^a): the values of the switching field are much higher, and more importantly, its behaviour as a 
function of the particle's size is opposite to that of the previous case, as now the field increases when the particle's 
size is decreased. For such high values of k s (K s » K c ) surface spins are aligned along their easy axes, and because 
of the strong exchange coupling they also drive core spins in their switching process. Thus the smaller the particle the 
larger the surface contribution, and the larger the field required for a complete reversal of the particle's magnetization. 
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FIG. 4. Switching field versus the surface anisotropy constant for ip = 0, j = 10 2 , and N — 10. 
In Fig. |] h c first slightly decreases for k s /j<0.l and then increases, and when k s becomes of the order of j it jumps 
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to large values. This plot illustrates the different regimes discussed above according to the value of the ratio k s /j. 
For an order of magnitude estimate of K s and the critical field H c at the point k s /j = 1, consider a 4 nm spherical 
cobalt particle for which J ~ 8 mev, K c ~ 2.7 x 10 6 erg/cm 3 . Then, K s ~ 10 erg/cm 2 , and H c ~ 5 T, when the field 
is applied along the core easy axis. 



III. CONCLUSION 



Considering the fact that experiments on nanoparticles show that the switching field does increase when the particle's 
size decreases (see e.g. || for cobalt particles), we may conclude that the anisotropy constant K s is at least of the 
order of the exchange coupling J, inasmuch as we can assume radial anisotropy on the surface, as is usually done in 
the literature. Then, as discussed above, for such values of K s , large deviations are observed with respect to the SW 
model in the hysteresis loop and thereby the limit-of-metastability curve, since in this case the magnetisation reverses 
its direction in a non-uniform manner via a progressive switching of spin clusters. So to deal with these features one 
has to resort to microscopic approaches such as the one used in this work. A work in progress applies the present 
technique to the cubo-octahedral cobalt particles with a diameter of about 3nm recently reported on in 0. 
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